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J> ' Abstract. We summarize the observations of the spin periods of rapidly accreting neu- 

tron stars. If gravitational radiation is responsible for balancing the accretion torque 
at the observed spin frequencies of w 300 Hz, then the brightest of these systems make 
excellent gravitational wave sources for LIGO-II and beyond. We review the recent the- 
' oretical progress on two mechanisms for gravitational wave emission: mass quadrupole 

. radiation from deformed neutron star crusts and current quadrupole radiation from 

' r-mode pulsations in neutron star cores. 



6 I SPINS OF ACCRETING NEUTRON STARS 



O ■ 

^ I Gravitational wave emission from rapidly rotating neutron stars (NS) has at- 

tracted considerable interest in the past several years. In addition to radiation 
from the spindown of newborn NSs (see the review by B. Owen in this volume), it 
rS I has long been suspected [1,2] that rapidly accreting NSs, such as Sco X-1, may be a 

I promising class of gravitational wave (GW) emitters. However, firm observational 

evidence of fast spins of these neutron stars had been missing until recently. 

NSs in low-mass X-ray binaries (LMXBs) have long been thought to be the pro- 
genitors of millisecond pulsars [3]. However, directly measuring their periods has 
proved elusive, probably because of their rather low magnetic fields. With the 
launch of the Rossi X-ray Timing Explorer, precision timing of accreting NSs has 
opened new threads of inquiry into the behavior and lives of these objects. RXTE 
observations [4] have finally provided conclusive evidence of millisecond spin pe- 
riods of NSs in about one-third of known Galactic LMXBs. These measurements 



^) Much theoretical progress has been made in understanding GW emission from LMXBs in the 
six months since the Amaldi conference in July 1999. Rather than just transcribe the talk given 
by one of us, we review the situation as of December 1999. Because of space limitations, this 
review is far from complete. 



8 



a 6 

D 
CO 

4 

O) 

a 

^ 9 







T 1 1 1 1 r 



Burst QPO 
kHz QPO I 
Pulsar I 




_l I I I I I I L 



200 400 600 800 
Frequency (Hz) 



10' 



10^ 



Frequency (Hz) 



FIGURE 1. Left (a) The distribution of spins of NSs in LMXBs, summarizing Tables 1-4 of [4]. 
Right (b) Characteristic signal amplitude he for several known LMXBs (symbols, see text) com- 
pared with the sensitivity /is/yr of LIGO-II in broadband (the dotted line) and narrowband (the 
solid line) configuration (provided by K. A. Strain on behalf of the LIGO Scientific Community). 
The dashed line shows LIGO-II sensitivity for a two- week integration (see text). 



are summarized in Fig. la. Altogether, there are seven such NSs with firmly estab- 
lished spin periods, by either pulsations in the persistent emission (discovered by 
Wijnands & van der Khs in the millisecond X-ray pulsar SAX J1808. 4-3658; [5]) or 
oscillations during type I X-ray bursts (burst QPOs, first discovered in 4U 1728-34 
by Strohmayer et al. [6]). There are an additional thirteen sources with twin kHz 
QPOs for which the spin may be approximately equal to the frequency difference 
[4]. A striking feature of all these neutron stars is that their spin frequencies lie 
within a narrow range, 260 Rz < Ug < 589 Hz. The frequency range might be 
even narrower if the burst QPOs seen in KS 1731-260, MXB 1743-29, and Aql X- 
1 are at the first harmonic of the spin frequency, as is the case with the 581 Hz 
burst oscillations in 4U 1636-536 [7]. These NSs accrete at diverse rates, from 
10-^^ Mq yr-i to the Eddington hmit, MEdd = 2 x lO^^M© yr'^ Since disk ac- 
cretion exerts a substantial torque on the NS and these systems are very old [8] , it 
is remarkable that their spin frequencies are so similar, and that none of them are 
near the breakup frequency of ^ 1.5 kHz. 

One possible explanation, proposed by White & Zhang [9], is that these stars 
have reached the magnetic spin equilibrium (where the spin frequency matches the 
Keplerian frequency at the magnetosphere) at nearly identical frequencies. This 
requires that the NS dipolar B field correlate very well with M [9,10]. However, 
there are no direct B field measurements for LMXBs, and in the strongly magnetic 
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binaries, where the B field has been measured directly, such a correlation is not 
observed. More importantly, for 19 out of 20 systems, there must be a way of 
hiding persistent pulses typically seen from magnetic accretors. These difficulties 
led Bildsten [11] to resurrect the conjecture originally due to Papaloizou & Pringle 
[12] and Wagoner [1] that gravitational radiation can balance the torque due to 
accretion. The detailed mechanisms will be discussed in the following sections. 

Regardless of the detailed mechanism for GW emission, if gravitational radiation 
balances the accretion torque, then it is easy to estimate the GW strength. As noted 
by Wagoner [1], in equilibrium the luminosities in GWs and in X-rays are both 
proportional to the mass accretion rate M, so the characteristic strain amplitude 
he depends on the X-ray fiux at Earth and the spin frequency 
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In Fig. lb we show he for Sco X-1 (marked with a star), a few other bright LMXBs 
with the spin inferred from kHz QPO separation (thick dots) and burst QPO fre- 
quency (triangles), and the miUisecond X-ray pulsar SAX J1808. 4-3658 (open di- 
amond). The dotted line shows LIGO-II sensitivity /i3/yr (i.e., h^ detectable with 
99% confidence in 10'' s, provided the frequency and the phase of the signal are 
known in advance [13]) in the broadband configuration, while the solid line shows 
/i3/yr for the narrowband configuration. However, the frequency and the phase are 
known precisely only for the SAX J1808. 4-3658 millisecond X-ray pulsar [14]. For 
other sources, Brady & Greighton [13] showed that the number of trials needed to 
guess the poorly known orbital parameters or to account for the torque noise due 
to M variations lowers the effective sensitivity by roughly a factor of two. 

While the average M certainly correlates with the X-ray brightness, current 
observations unfortunately do not let us robustly infer the instantaneous torque [4] . 
Even though M varies on a timescale of days, torque noise leads to frequency drift 
only on a timescale of weeks. The accretion torque is A^a — M{GMRY/'^ , and the 
total time-averaged torque is zero due to equilibrium with GW emission. Assume 
that Na flips sign randomly on a timescale tg ~ few days. The spin frequency Vt 
will experience a random walk with step size 8Vl = {Na/I)ts, where I is the NS 
moment of inertia. After an observation time tobs) the drift is AQ = (tobs/^s)^^^^^- 
This will exceed a Fourier frequency bin width, i.e., AQ ^ 27r/tobs only after 

t^^^ ^ 21 days / Iday \ / iQ-' yr-' \ ^'^ 



Hence, on a timescale of tens of days, the intrinsic GW signal is coherent. The 
dashed line in Fig. lb shows the LlGO-11 sensitivity for a two-week integration in a 
narrowband conflguration. This suggests that the way to detect GWs from LMXBs 
may be short integrations [13]. 

Gurrently, there are two classes of theories for GW emission from NSs in LMXBs. 
The presence of a large-scale temperature asymmetry in the deep crust will cause 
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it to deform [11]. The resulting "mountains" will give the rotating star a time- 
dependent mass quadrupole moment. Alternatively, unstable r-mode pulsations 

(see a review by B. Owen in this volume) of a suitable amplitude in the NS liquid 
core can emit enough gravitational radiation to balance the accretion torque [11,15]. 

II DEFORMATIONS OF ACCRETING NS CRUSTS 

The crust is a 1 km layer of crystaUine "ordinary" (albeit neutron-rich) matter 
that overlies the liquid core composed of free neutrons, protons, and electrons. The 
crust's composition varies with depth in a rather abrupt manner. As an accreted 
nucleus gets buried under an increasingly thick layer of more recently accreted ma- 
terial, it undergoes a series of e~ captures, neutron emissions, and pycnonuclear 
reactions [16-18], resulting in layered composition. In Fig. 2a, we show schemati- 
cally two such compositional layers (light and dark shading) sandwiched between 
the liquid core and the ocean. Since an appreciable fraction of the pressure is sup- 
plied by degenerate electrons, e~ captures induce abrupt density increases. In the 
outer crust, these density jumps are as large as ~ 10%, while in the inner crust 
the density contrast is ^ 1%. At T = 0, the e~ captures occur when the elec- 
tron Fermi energy E-p is greater than the mass difference between the e~ capturer 
and the product of the reaction. In the absence of other effects, this depth is the 
same everywhere, and such an axisymmetric capture boundary (the dashed line in 
Fig. 2a) does not create a mass quadrupole moment. 

However, in accreting NSs the crustal temperatures are high enough (in excess 
of 2 X 10^ K) that e^ capture rates become temperature-sensitive [19]. Bildsten 
[11] pointed out that if there is a lateral temperature gradient in the crust (the 
arrow in Fig. 2a), then regions of the crust that are hotter undergo captures at 
a lower density than the colder regions. The capture boundary becomes "wavy" 
(the solid line in Fig. 2a), with captures proceeding a height Az^ higher on the hot 
side of the star, and A^;^ lower on the cold side. Such a temperature gradient, if 
misaligned from the spin axis, will give rise to a nonaxisymmetric density variation 
and a nonzero quadrupole moment Q22 [H]- 

The required quadrupole moment Qeq such that GW emission is in equilibrium 
with the accretion torque is 

The range of M's in LMXBs is f» 10~^^ — 2 x 10"^ Mq yr"\ requiring Q22 ~ 
10^^ — 10^^ g cm^ for i/s = 300 Hz [11]. Can temperature-sensitive e~ captures 
sustain a quadrupole moment this large? The quadrupole moment generated by 
a temperature-sensitive capture boundary is Q22 ~ Qm = ApAzj-R^, where Ap is 
the density jump at the electron capture interface. Qgd is the quadrupole moment 
that would result if the crust did not elastically adjust (or just moved horizontally) 
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FIGURE 2. Left (a) A cartoon description of how a transverse temperature gradient leads to a 
varying altitude for the e~ captures. Right (b) The quadrupole Q22 due to a single capture layer 
in the inner crust as a function of Af. The solid and dashed lines denote the results for several NS 
models for a 10% composition asymmetry, while the dotted line is the relation given by Eq. (3), 
i.e., the quadrupole necessary for spin equilibrium at Vs = 300 Hz as a function of M. 



in response to the lateral pressure gradient due to wavy e~ captures. Using this 
estimate, Bildsten [11] argued that a single wavy capture boundary in the thin outer 
crust could generate Q22 sufficient to buffer the spinup due to accretion (Eq. [3]), 
provided that temperature variations of ~ 20% are present in the crust. 

However an important piece of physics is missing from this estimate: the shear 
modulus /i. If yU = 0, the crust becomes a liquid and cannot support a non-zero 
Q22- Ushomirsky, Cutler, & Bildsten [20] recently calculated of the elastic response 
of the crust to the wavy e~ captures. They found that the predominant response 
of the crust to a lateral density perturbation is to sink, rather than move sideways. 
For this reason, Q22 generated in the outer crust [11] is much too small to buffer the 
accretion torque. However, a single e~ capture boundary in the deep inner crust 
can easily generate an adequate Q22- Because of the much larger mass involved by 
captures in the inner crust, the temperature contrasts required are <^ 5%, or only 
fa 10^ — 10'^ K, not ^ 10^ K as originally postulated [11]. 

What causes the lateral temperature asymmetry, and can it persist despite the 
strong thermal contact with the almost perfectly conducting core? In LMXBs, the 
crusts are composed of the compressed products of nuclear burning of the accreted 
material. The exact composition depends on the local accretion rate, which could 
have a significant non-axisymmetric piece due to, e.g., the presence of a weak B 
field. Moreover, except in the highest accretion rate LMXBs, nearly all of the 
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nuclear burning occurs in type I X-ray bursts. Burst QPOs (see Sec. I) provide 
conclusive evidence that bursts themselves are not axisymmetric. Until the origin 

of this symmetry breaking is clearly understood, it is plausible to postulate that 
these burst asymmetries get imprinted into the crustal composition. 

Ushomirsky et al. [20] showed that such a non-uniform composition leads directly 
to lateral temperature variations 5T. Horizontal variations in the charge-to-mass 
ratio Z"^ I A (which determines the crustal conductivity) and/or nuclear energy re- 
lease modulate the radial heat flux in the crust and set up a nonaxisymmetric bT . 
The 5T's required to induce a Q21 ~ Qeq can easily be maintained if there is a 

10% asymmetry in the nuclear heating or Z"^ jA. So long as accretion continues, 
these (5T's persist despite the strong thermal contact with the isothermal NS core. 

The e~ capture Q22 calculations [20] are summarized in Fig. 2b. If the size of 
the heating or Z"^ j A asymmetry is a constant fixed fraction, then for M <^ O.SMsdd 
the scaling of Q22{M) is just that needed for all of these NSs to have the same spin 
frequency (the normalization is proportional to the magnitude of the asymmetry, 
but the scaling is fixed by the microphysics) . For M ^ O.SMEdd; in order to 
explain the spin clustering at exactly 300 Hz, this mechanism requires that the 
crustal asymmetry correlate with M. Alternatively, if the asymmetry is the same 
as in the low M systems, then one would expect the bright LMXBs to have higher 
spins, a possibility that cannot be ruled out by current observations (Sec. I). 

So long as crustal deformations are due to shear forces only, the crustal Q22 is 
limited by the yield strain a„iax to be less than [20] 

\ p6.26 

(522's needed to buffer the accretion torque require strains a ~ 10"'^ — 10~^ at 
300 Hz, with a ^ 10~^ in near-Eddington accretors. Estimates for the yield 
strain of the neutron star crust range anywhere from 10~^ for perfect one-component 
crystals to 10^^. Hence a ^ 10^^ is probably higher than yield strain, though this 
conclusion is based on extrapolating experimental results for terrestrial materials 
by > 10 orders of magnitude. Such high strains are perhaps the biggest problem 
with the crustal Q22 mechanism. At high pressures (^ shear modulus) terrestrial 
materials tend to deform plastically rather than crack, and so the crusts of accreting 
NSs may be in a state of continual plastic flow. If accretion continually drives the 
crust to cXmax, this leads to a natural explanation for spin similarities near MEdd- 

However, many fundamental issues remain unanswered. First, the calculation 
[20] is only good up to an overall prefactor set by the density of capture layers in 
the deep crust. We thus need an exploratory calculation of both the composition 
of the products of nuclear burning in the upper atmosphere over the entire range of 
M in LMXBs, and their detailed nuclear evolution under compression in the crust. 
Knowledge of the composition is also necessary for a robust calculation of the shear 
modulus, which is clearly the crucial number to know when computing the elastic 
response of the crust. Recent results [21,22] indicate that inner crusts of NSs are 
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composed of highly nonspherical nuclei and may be more like liquid crystals (solids 
that provide no elastic restoring force for certain kinds of distortions) rather than 

simple Coulomb solids. Such improved calculations have implications far beyond 
the problem of the crustal quadrupole moment. The shear modulus of the crust 
affects the maximum elastic energy that can be stored in the crust, and hence the 
energetics of pulsar ghtches and starquakes, as well as the models of magnetic field 
evolution that depend on crustal "plate tectonics" (see [23]). It even has bearing on 
the stability of r-modes in neutron stars (Sec. 111). In addition, much work needs 
to be done on understanding what sets the shear strength cTmax of multicomponent 
crystals, likely with defects and highly nonspherical nuclei, or what happens when 
^max is exceeded and viscoelastic flow ensues. 

Ill R-MODES IN ACCRETING NS CORES 

Bildsten [11] and Andersson, Kokkotas, & Stergioulas [15] pointed out that the 
r-mode instability (see the review by B. Owen in this volume for the introduction 
and notation) may also explain the spins of NSs in LMXBs, and, if so, produce 
GW signal detectable by LIGO-II. An accreting NS is spun up (along a line in 
(i^s, T) plane marked with an arrow in Fig. 3) until it reaches the r-mode instability 
line (the solid line in Fig. 3). At that point (marked by a thick dot in Fig. 3) 
the r-mode amplitude needed to balance the accretion torque is rather small. The 
NS can then hover at the instabihty hue, with I/tg + 1/ry = 0, and the r-mode 
amphtude such that it balances the accretion torque. However, at T = few x 10^ K, 
the r-mode— accretion equilibrium spin frequency would he ^ 150 Hz, rather than 
~ 300 Hz, resulting in an apparent disagreement with the observed spins of LMXBs. 
Bildsten [11] and Andersson et al. [15] speculated that including other sources of 
viscosity, e.g., superfluid mutual friction, is likely to raise the instability curve, 
resulting in equilibrium frequencies closer to the canonical 300 Hz. Finally, the 
narrow range of the observed spin frequencies would presumably arise because of 
the similar core temperatures of the accreting NSs (shown by the shaded box in 
Fig. 3). 

Recent theoretical work brought up several challenges to this scenario. Levin 
[24] and Spruit [25] showed that steady-state equilibrium between accretion and 
r-modes is thermally unstable for normal fluid cores. In a normal fluid (i.e., not 
superfluid), the shear viscosity scales as T^^, so the increase in the core temperature 
due to viscous heating decreases the shear viscosity. The smaller shear viscosity 
increases the growth rate of the r-mode, leading to an unstable runaway. Using 
a phenomenological model of nonlinear r-mode evolution [26], Levin [24] showed 
that in this case, instead of just hovering near the instability line, the r-mode 
grows rapidly until saturation, heats up the star, and spins it down and out of the 
instability region in less than 1 yr. Therefore, if NSs in LMXBs have normal fluid 
cores, we would not expect to see any of them with ^ 300 Hz spins. 

The unstable regime for r-modes in normal fluid NSs (above the solid line in 

7 



iiUl I I I 

10^ 10^ 10^ 10^" 10^^ 
Core Temperature (K) 

FIGURE 3. Critical spin frequencies of the r-mode instability. The solid line is the critical 
frequency set only by shear and bulk viscosities in the core (same as Fig. 1 of B. Owen's review in 
this volume). The dashed line also includes viscous boundary layer damping, while the shading 
around it displays the effect of core superfluidity. 

Fig. 3) encompasses much of the parameter space occupied by NSs in LMXBs and 
newborn NSs. Because of the large torques exerted by the unstable r-modes, we 
would not expect to see any NSs in this region. In addition, the existence of two 
1.6 ms radio pulsars (the spins and upper limits on core temperatures of which are 
shown by arrows in Fig. 3) means that rapidly rotating NSs are formed in spite 
of the r-mode instability [15]. While it is not clear whether their current core 
temperatures place these pulsars within the r-mode instability region, normal-fluid 
r-mode theory says that they were certainly unstable during spinup. 

Superfluid r-mode calculations have been eagerly awaited, as they could resolve 
these conflicts. However, Lindblom & Mendell [27] showed that, for most val- 
ues of the neutron-proton entrainment parameter, the superfluid dissipation is not 
competitive with gravitational radiation. Only over about 3% of the possible en- 
trainment parameter values is mutual friction strong enough to compete with grav- 
itational radiation. The r-mode instability line in this case is an approximately hor- 
izontal line (see Fig. 8 of [27]) separating the unstable spin frequencies {ug > i^cmt) 
from the stable ones {ug < i^ait)- If the superfluid entrainment parameter has a 
value such that z/crit ~ 300 Hz, then the LMXB spin frequencies could still be 
understood in terms of the r-mode instability and the special nature of the NS 
superfluid. 

Before learning about these results. Brown & Ushomirsky [28] ruled out such a 
simple superfluid equilibrium observationally for a subset of LMXBs. In steady 



state, the shear in the r-mode deposits f« 10 MeV of heat per accreted baryon 
into the NS core. When the core is superfluid, Urea neutrino emission from it is 

suppressed, and this heat must flow to the NS surface and be radiated thermally. 
In steadily accreting systems (such as Sco X-1) this thermal emission is dwarfed by 
the accretion luminosity of GMM/R ^ 200 MeV per accreted baryon. However, in 
transiently accreting systems, such as Aql X-1, when accretion ceases, the r-mode 
heating should be directly detectable as enhanced X-ray luminosity from the NS 
surface. For Aql X-1 and other NS transients. Brown & Ushomirsky [28] showed 
that, if the superfluid r-mode equilibrium prevails, then the quiescent luminosity 
should be about 5—10 times greater than is actually observed. 

A possible resolution of this conundrum has been recently proposed by Bildsten 
& Ushomirsky [29]. All but the hottest 10^° K) NSs have sohd crusts. The 
r-mode motions are mostly transverse, and reach their maximum amplitude near 
the crust-core boundary. The fluid therefore rubs against the crust, which creates 
a thin (few cm) boundary layer. Because of the short length scale, the dissipation 
in this boundary layer is very large. The damping time due to rubbing is [29] 

Trub ~ 100 S , (5) 

^6 \ / 

substantially shorter than the viscous damping times due to the shear and bulk 
viscosities in the stellar interior, as well as the mutual friction damping time for 
most values of the superfluid entrainment parameter. 

The critical frequency for the r-mode instability in NSs with crusts is shown by 
the dashed line in Fig. 3 for the case where all nucleons are normal, and the dark 
shading around it represents the range of frequencies when either neutrons or all 
nucleons are superfluid. The crust-core rubbing raises the minimum frequency for 
the r-mode instabihty in NSs with crusts to ^ 500 Hz for T ^ 10^° K, nearly 
a factor of five higher than previous estimates. This substantially reduces the 
parameter space for the instability to operate, especially for older, colder NSs, such 
as those accreting in binaries and millisecond pulsars. In particular, the smallest 
unstable frequency for the temperatures characteristic of LMXBs is ^ 700 Hz, safely 
above all measured spin frequencies. This work resolves the discrepancy between 
the theoretical understanding of the r-mode instability and the observations of 
millisecond pulsars and LMXBs, and, along with observational inferences [28], likely 
rules out r-modes as the explanation for the clustering of spin frequencies of neutron 
stars in LMXBs around 300 Hz. 

To summarize, a significant role of steady-state r-modes in LMXBs has proba- 
bly been ruled out, both on theoretical grounds [29] (unless crust-core coupling is 
much stronger than was estimated), and observationally [28] (for Aql X-1 in par- 
ticular). However, stochastically excited r-modes that decay rapidly may still play 
a significant role in accreting systems, as even a very small amplitude (a ^ 10~^) 
can balance the accretion torque at MEdd- In addition, the issues of crustal shear 
modulus and the structure of the crust-core boundary, highlighted in Sec. II, are 
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of paramount importance for r- modes as well. Crustal quadrupoles [11,20] can ex- 
plain the spins of LMXBs and remain a viable source of continuous GWs, but the 

strains at MEdd are rather high. The crustal breaking strain cTmax is not likely to be 
understood theoretically any time soon, and detection of GWs from LMXBs with 
LIGO-11 type instruments will surely teach us many new things about NSs. 
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